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Abstract We study the velocity distribution of unicellular swimming algae Eu¬ 
glena gracilis using optical microscopy and theory. To characterize a peculiar fea¬ 
ture of the experimentally observed distribution at small velocities we use the 
concept of active fluctuations, which was recently proposed for the description 
of stochastically self-propelled particles [Romanczuk, P. and Schimansky-Geier, 
L., Phys. Rev. Lett. 106 , 230601 (2011)]. In this concept, the fluctuating forces 
arise due to internal random performance of the propulsive motor. The fluctuat¬ 
ing forces are directed in parallel to the heading direction, in which the propulsion 
acts. In the theory, we introduce the active motion via the depot model [Schweitzer 
et al, Phys. Rev. Lett. 80 , 23, 5044 (1998)]. We demonstrate that the theoretical 
predictions based on the depot model with active fluctuations are consistent with 
the experimentally observed velocity distributions. In addition to the model with 
additive active noise, we obtain theoretical results for a constant propulsion with 
multiplicative noise. 
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1 Introduction 

The development of the theory of Brownian motion, connected to such famous 
names as Einstein [1], Langevin [2] or Smoluchowski [3], was one of the greatest 
triumphs in statistical physics at the beginning of the 20*** century. Since then 
Brownian motion not only assumed a central role in the foundation of thermody¬ 
namics and statistical physics but continues to date to be a major interdisciplinary 
research topic. 

Recently there has been an increasing focus in biology, physics and chemistry, 
on so-called active matter systems, both from experimental and theoretical point 
of view. Examples of such far-from-equilibrium systems range from the dynami¬ 
cal behavior of individual units such as Brownian motors [1], motile cells [SHini, 
macroscopic animals [nHi3] or artificial self-propelled particles [niGHi to large 
coupled ensembles of such units and their large scale collective dynamics |19H26| . 
Major advances in active matter research have been made possible by the contin¬ 
uously improving experimental techniques such as for example automated digital 
tracking H El na [23] or the realization of active granular and colloidal systems 
[Il[l5l|Il|27]. 

Despite these recent advances, there are still many open questions regarding 
the universal properties of active motion. Already at the level of individual ac¬ 
tive particles it is important to distinguish thermal fluctuations resulting from 
the interactions with a heat bath from effective stochastic forces which have their 
origin in the active nature of the system. Recently, it was shown that such active 
fluctuations may lead to a complex behavior of the mean squared displacement of 
self-propelled particles with multiple crossovers esehi, and may manifest them¬ 
selves in characteristic behavior of probability density functions at small speeds, 
strikingly different from the ones obtained for thermal fluctuations |251 [29] . Non¬ 
trivial motion statistics and patterns also arise in systems that can switch between 
two or more locomotion modes like run-and-tumble [sniEi] or two different speed 
states EiiEa. 

We should note that the term “active fluctuations” has been used before to 
refer to non-equilibrium fluctuations of intracellular filament networks (see e.g. 
[33H35] f. and recently also in the context of flagellar beating [36]. Although all 
these concepts can be related on a very fundamental level, here we restrict our 
definition of active fluctuations to the velocity fluctuations of self-propelled agents, 
which can be described as active Brownian particles |25| . 

In this work we report on experimental observed velocity distribution of uni¬ 
cellular swimming algae Euglena gracilis, which cannot be explained by models 
with passive noise only. To describe the motion of the cells, we develop a theory 
of active Brownian motion with active fluctuations and derive the corresponding 
particle velocity distributions. We demonstrate that the measured velocities are 
consistent with the theoretical predictions, which allows for insight into the motion 
characteristics and control patterns. 

The remainder of the paper is organized as follows: In Section we describe 
the experimental setup used to measure velocity distributions of Euglena gracilis. 
Further on, in Section]^ we introduce the theory and discuss the comparison be¬ 
tween the theory and experiment. In addition, we discuss also a second theoretical 
model, to demonstrate that multiplicative noise can be also considered within the 
active fluctuations framework. Finally, we conclude in Section]^ 
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2 Experiment 

2.1 System 

In this work, we studied the motion of unicellular photosynthetic flagellate Eu¬ 
glena gracilis Klebs using optical microscopy. It is known that Euglena gracilis is 
most sensitive to blue light: the photosynthetic efficiency is maximal for 465 nm 
light m- Euglena is a flagellated swimmer (a pusher-type), propelled by a single 
flagellum. The motion is not purely phototactic: Euglena cells in dark conditions 
perform random motion with much spinning |38j . The threshold intensity of blue 
light required to cause the negative phototaxis is approximately 10 mW/cm^ [39] . 
This is an almost two times lower limiting light intensity than that for a normal 
daylight m- In our experiments we investigated behavior of the cells in strong 
blue light beyond the photosynthetic limit when the photophobic behavioural reac¬ 
tions are most pronounced. The cells were supplied by Carolina Biological Supply 
Company (Burlington, NC, USA) and grown on an aerated, inorganic salt medium 
at room temperature under a 12:12 hours light-dark cycle. The average cell size 
(length) was about 40 pm. 

The schematic of the experimental setup is shown in Fig.[^ During each exper¬ 
imental trial the cells were confined inside a microfluidic channel (50 mm length, 
5 mm width, 400 pm height and 100 pi volume) of an ibiTreat 0.4 Luer channel 
slide installed into a slide-mount on the stage of the inverted microscope. Small 
height of the channel allowed us to consider motion in two spatial dimensions. The 
tubes on both ends of the channel were kept open during the experiment ensuring 
proper gaseous cell metabolism. Illumination of the cells with light was performed 
by two light emitting diodes (LED) positioned in the central part of the chan¬ 
nel so that the angle between the longitudinal axis of each LED and the channel 
plane was fixed at 60° and the distance between the centers of diodes was 10 mm. 
This arrangement takes into account characteristics of the directivity pattern of 
each LED and ensures uniformity of the illumination within the view field of the 
optical system. In this work, we used a pair of blue LEDs with a peak emission 
at 465 nm (NSPB500S, Nichia Corporation , Japan). The diodes were lighted up 
simultaneously. The light intensity was set to 50 mW/cm^ - well beyond the pho- 
tosynthetic limit. The activation light from the diodes has been used for detection, 
no additional illumination was required. All experiments were performed at room 
temperature. Duration of each run was about 10 minutes. 


2.2 Imaging and data collection 

The optical readout system consisted of an image acquisition equipment and data 
processing computer. Real time images of Euglena cells in a channel slide were 
captured by the Andor iXonS 885 EMCCD camera (1004 x 1004 pixels, pixel 
size 8 pm, 2 fps frame rate) mounted on the optical microscope Zeiss Axiovert 
200M through a 5x objective lens. The observation field was 1.6 mm x 1.6 mm. 
The number of cells within the observation area was 800 — 1200 giving packing 
fractions of 0.026 — 0.039. Images were recorded in miltipage TIFF format using 
original camera manufacturer software. Typical trajectories of the cells are shown 
in Fig. 
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Fig. 1 Scheme of the experimental setup. 



Fig. 2 Typical trajectories of Euglena gracilis cells captured in strong blue light. Only ten 
trajectories are shown for the sake of clarity. 


Further image processing involving cell recognition and tracking was performed 
with Imaris package, v7.6 (Bitplane, UK). The raw data consisted of x and y coor¬ 
dinates of the cell geometrical centers, identities of cells and a time stamp. These 
data did not include positions of non-motile (dead) cells which were excluded with 
a filter based on a minimal trajectory length (50 pm). In the videos, the frame 
rate seems to be sufficiently high. The displacements between the frames are less 
than the cell size. The positional statistics was subsequently exported into ASCII 
format files and analyzed with our code. Velocities of the cells were calculated 
based on the two subsequent positions. The motion statistics were averaged over 
two runs. Due to isotropy of the motion, we combined the histograms in x and y 
directions. To extract the slice of the histogram corresponding to one velocity com¬ 
ponent {wx,0), we used the data where the corresponding transverse component 
Vy was small —0.5 pm/s < Vy < 0.5 pm/s. 

The cells in strong blue light were in constant motion, although slowing down 
from time to time. As we noticed the change of the motion within the duration 
of the experiment, we split the statistics into 6 subsets corresponding to time 
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Fig. 3 (a) An example of raw data set for Euglena gracilis in strong blue light collected 200- 
300 seconds from the start of the illumination, (b) Six symmetrised datasets, where and 
■Evy are combined. 


intervals of 0-100, 100-200, 200-300, 300-400, 400-500, and 500-600 seconds. An 
example of a raw dataset with Vx and Vy velocity components for one interval and 
the collection of symmetrised data are shown in Fig. We can observe straight 
away a striking feature of these distributions: A strong peak around zero velocity. 
This is in contrast with the shape of distribution predicted for active Brownian 
particles effected by thermal agitation, likewise usual Brownian motion which have 
a crater like shape with a minimum around zero. This agitation from the molecules 
of the surrounding liquid is negligible compared to the phototactic propulsion. To 
understand the peak at vanishing velocity we need to introduce a different noise 
model, which we present in the next section. 


3 Self propelled particles with active fluctuations 

3.1 Models with active fluctuations 

In recent publications two of us have investigated the difference between passive 
and active noise acting on the motion of self propelled units |251I29| . Passive fluctu¬ 
ations have their origin in the fluctuating environment in which the particle moves. 
In result, the effect of a passive random force » 7 p(t) is independent on the direction 
of motion (heading) of the particle. The classical example of passive fluctuations 
is ordinary Brownian motion in thermal equilibrium, where the stochastic force is 
associated with random collisions of the Brownian particle with molecules of the 
surrounding fluid. 

In contrast, active fluctuations are a pure far from equilibrium phe¬ 

nomenon. They may become important in the motion of an biological agents 
or self-propelled particles. The origin of these fluctuations can be for example 
variations in the propulsion of chemically powered colloids [me], complex intra¬ 
cellular processes in cell motility 13 [8] or unresolved internal decision processes in 
animals |11H13| . 

In this manuscript we show that the concept of active fluctuations is able to 
explain measurements of the velocity distribution function of the Euglena cells. 
Especially, an observed giant peak in the two dimensional velocity distribution of 
the Cartesian velocity components P{vx,Vy) at the origin Vy = Q,Vy = Q is due to 
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the special selection of the noise sources. Only active fluctuations describe correctly 
the large peak in two-dimensional Cartesian velocity distributions at small speeds. 

We assume a two-dimensional motion of the cells is given by the position vector 
r{t) and the velocity v{t). Components are labeled as r = {x,y) for the position 
vector of the cell, and v = {vx,Vy) being its velocity vector. The dynamics of the 
kinematics of the cell is given by the second Newton law. We assume a propulsive 
mechanism with force ap, Stokes friction and a random forces ri{t) acting on the 
cell. In particular, the dynamics reads afterwards. 

r = V 

mv = ap — 7 V -|- ri{t) 

and we will set the mass to unity m = 1, later on. 

Active particles have a preferred direction of motion (heading) determined by 
their propulsion mechanism (“head-tail” asymmetry). The heading vector is a unit 
vector Bfi with |e/i| = 1 dehning the orientation of the particle along which the 
propulsive force acts. This force is independently on other acting external forces. 
It acts generally in direction of the instantaneous velocity, and can be decelerating 
or accelerating. 

In two spatial dimensions {d = 2) it is fully determined by the angle defining 
the direction with respect to the s-axis (see Fig. |4|i): e/i(t) = {cosifi{t),siinp{t)), 
and the evolution of the position of the particle can be rewritten as 

r{t) =v{t)eh{t). (3) 

Note, that v{t), which is the projection of the velocity onto the heading direction, 
can also adopt negative values if moving backwards with respect to the propulsive 
mechanism. The heading vector differs from the unit vector of the velocity (t) to 
which it is either parallel or anti-parallel if moving forward or backward. Projected 
on 01 , (t), the velocity is the speed and always positive and the length of the vector 
v(t) = s{t)ev{t) with s{t) = \/vx + Vy = |u(t)|. 

The temporal evolution of the velocity vector in the new coordinates (u,(p) 
reads 


( 1 ) 

( 2 ) 


V = veh+ vipe^ , (4) 

where = {—sin cos is the unit vector in the angular direction. Multi¬ 
plying Eq. (2) with e„ and e<^, respectively, yields: 


v = ap-'^v + 'n-eh, 

(5) 

1 

ip = -T) ■ Bip . 

(6) 


Here, • indicates a scalar product. One has to note that the orientational dynamics 
diverge for u = 0. This is due to the fact that for a point-like particle, as considered 
here, the polar angle (p cannot be dehned for a vanishing speed. 

At this point one has to distinguish between the action of active vs passive 
noise. Passive noise as in case of Brownian motion acts independently in the x and 
y directions ex, respectively Sy, as Gaussian white noise. It decomposes into two 
independent components with 

T/p(t) = y/W[^x{t)ex + iv{t)ey), 


(7) 
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Fig. 4 (a) Schematic visualization of motion in a plane with unit vectors e^(t),e<p(t) (thick 
arrows), velocity vector v(t) (thin arrow) and dsh/dl (thick red/gray arrow). The trajectory of 
the particle l{t) is indicated by the dotted line, (b) Visualization of passive fluctuations ix{t) 
and ^y{t) (thick red/gray arrows), (c) Visualization of active fluctuation ^v{t) and (thick 

red/gray arrows). Figure reproduced from [49]. 


with following properties 

= , i,je{x,y}. (8) 

Thus with this choice of velocity variables, the problem becomes a model with 
multiplicative noise in D- It is easier solved in a Cartesian representation of the 
velocity vectors, where for passive noise the motion in x and y becomes indepen¬ 
dent. 

In contrast to passive noise, for active fluctuations we assume Gaussian noise 
acting with intensity Dy along the heading direction. It is complemented by a 
second Gaussian noise randomly changing the heading direction with intensity 
D^p, by acting only perpendicularly to the heading vector. Therefore, we define 

'n^{t) = V2Dv^v{t)eh{t) + (9) 

Again both noise sources shall be independent 

(Ci(^))=0. , i,je{v,ip}. (10) 

A similar decomposition happens in case of active noise compared to the Carte¬ 
sian formulation with passive noise. The dynamics along the heading direction 
with velocity component v decouple from the angular motion. This stochastic </?- 
dynamics with homogeneous additive noise along circles of constant v results even 
in an equipartition along ip. Subsequently the v dynamics can be solved. As a re¬ 
sult, models with active fluctuations becomes a new class of analytically tractable 
systems. 
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3.2 Energy depot 


To explain measnrements of the velocity distribution function of the Euglena cells, 
we will use the depot model as introduced by Ebeling, Schweitzer and Tilch [251 mi¬ 
lls] . This model traces the propulsion back to a food uptake of the cells which rate 
q from the environment. It fills up an internal energy depot Cinif) of the cell. The 
internal energy decay with rate c and can be transformed into kinetic energy with 
rate d, whereby the corresponding term is proportional to the internal and kinetic 
energy of the cell. Thus ei„(t) obeys 

Cm — Q dc^nV . (f f) 


The gained kinetic energy is transformed into momentum and creates a propul¬ 
sion in the Newtonian force balance along the heading direction as 

V = de{t)v{t) — + V2JX^v{t). (12) 


We assume that the angular dynamics are purely stochastic, driven by Gaussian 
white noise: 

(13) 

This assumption implies also that the angular dynamics are not affected by the 
depot. In general, this might be not the case, for example angular fluctuation could 
decrease with decreasing energy depot. However, this would only affect the per¬ 
sistence of the cell trajectory and would make no difference to stationary velocity 
distributions we focus on here. The active noise sources ^fe(t) are defined in Eq. 

Assuming a fast relaxation rate of the depot variable, the internal energy re¬ 
sponds instantaneously to its supply and demand. Thus it assumes a stationary 
value and can be adiabatically eliminated, despite ongoing exchange with the envi¬ 
ronment (for a further discussion including stochastic supply q, see jS^). It yields 


e(t) = 


c + dv‘^{t) ’ 

and in consequence the heading velocity changes as 


dq 


- _ 'yu -{- y/O.Dv^v (t) . 


c -h dv'^ 

Hence, the propulsive acceleration reads in this approximation 

dq 

^ c -h 

Einally, the depot model becomes 

V = CLp 'yv (f)! ^ ~ ^ (t). 


(14) 


(15) 


(16) 


(17) 


and the corresponding Eokker-Planck equation 


—Pa(v,¥3,t|i;o,<^o,to) 


d_ 

dv 


(flp — 


7«)pa + Dv 


dv^ 


Pa + 


9 r\ 9 • 


(18) 
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with Op from Eq. (16). 

The propulsion dynamics can be reformulated by introducing the stationary 
propulsion speed in the absence of noise vq defined as 


q c 


Eq. (15) becomes then 


'yd 


){vo - v^) + y/Wviv{t). 


(19) 


( 20 ) 


c + dv'^ 

In case that parameters obey ® ^ > 0, uq exists and is a stable fixed point of 

the deterministic speed dynamics. Its positive root wq > 0 is the mean propulsion 
speed in the small noise limit. The depot pumps sufficient energy to excite a 
stationary motion. In the opposite case, the single stationary stable solution of the 


deterministic part of Eq. (15) is the rest state (vg = 0). The depot is insufficient 
to generate persistent motion at the given damping 7 . Thus, in this case particles 
are driven by the noise alone. Note that this noise has to be still associated with 
active fluctuations produced by the propulsion mechanism of the algae as it acts 
only along the heading direction, which is fundamentally different from the case 
of thermal agitation. 


Later on, we estimate the value of vg from Eq. (19) using experimentally fitted 
parameters. 


3.3 Additive active noise 


The asymptotic stationary distribution limtg->.-ooPa(v,(p,tlvg,tpg,tg) = pa(v) be¬ 
comes independent of the initial conditions. Also the noisy v(t) dynamics become 
independent from the direction ip and in the stationary limit the probability func¬ 
tion density decomposes pa(v,ip) = pa(v)Pd(‘P)- The angular distribution becomes 
homogeneous in all directions p with Pd{p>) = l/27r. The remaining probability 
density function for the velocity follows from the stationary Eokker-Planck equa¬ 
tion 


dt 


Pa{v) = 0 


d_ 

dv 


(ap — 


yvjpa + Dv 




( 21 ) 


with Op from Eq.(16). The latter ordinary differential equation can be integrated 
with vanishing probability flux. Thus, the stationary mobability density function 
of velocities projected on the heading dynamics reads [j 


Pa 




( 22 ) 


Erom this, by summing the positive and negative projected velocities, we obtain 
the speed distribution, i.e. p(s) = pa(s = v) + pa(s = —v) It reads 


p(s) =Pa(s 


v) +Pa(s 





(23) 


^ Probability density functions are defined up to the normalization constant. 

2 Note, it follows from this definition that the derivative of the speed distribution p(s) 
vanishes in the limit s —> 0 if pa is differentiable at u = 0. 
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Fig. 5 T he stationary speed probability (Eg. | |23^ , left panel) and velocity distribution function 
(Eg. ( |24[ |, right panel) for the energy depot with additive noise. Parametersig = 2,d = l,c = 
1,7 = 1 and noise intensity as assigned in the inset. 


Transformation into the Cartesian velocities requires the consideration of the 
Jacobi-determinant for the two-dimensional variable transformation s,i/p —^ Vx,Vy. 
Therefore a giant peak occurs at the origin reflecting the non-vanishing speed prob¬ 
ability density for the “resting” situation in the presence of active noise. Finally, 
we obtain 


P{Vx,Vy) oc 




vl + vl 


/ c 2 , 2V 


exp 



(24) 


The divergence of the Cartesian probability density is a universal property of 
active noise acting along the heading direction [251129] . It is intrinsically connected 
to finite values of the speed probability density at the origin p{s = 0) > 0. This in 
turn, is a consequence of active noise stochastically exploring speed states infinitely 
close to the resting situation. These states are never reached in case of an ordinary 
Brownian particle, where the permanent stochastic impacts, independent in x- and 
^-direction, always drive the speed dynamics away from the state s = 0 . 


The solution Eq. \22\ looks very similar to the stationary solution in the depot 
model as reported in [44] . However it has to be emphasized that it describes a 
different physical situation. The solution in [13] describes the Cartesian velocity 
components as the results from a passive noise as given by Eq. 0. In the case 
discussed here the solution was obtained along the co-moving heading direction 
and Cartesian velocity components are due to Eq. (24). 


As shown in a previous paper m the existence of the state with a maximum 
in the probability density function can act as a strong attractor in a situation 
with many particles and a velocity-alignment interaction between particles. In this 
situation the particles with active noise are collectively “trapped” in the resting 
state in contrast if passive noise is applied. The character of the phase transition 
between the ordered (running particles) state and disordered state changes from a 
second order to a first order transition due to the coexistence of the running and 
resting states. 
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Fig. 6 Euglena gracilis one-di men sional velocity histograms collected in intense blue light 
(points) and the fits with Eq. ( |24[ l (lines). The labels on the experimental sets indicate the 
time interval within which the data were collected: 100 corresponds to the 0-100 seconds, 300 
to 200-300 seconds, and 500 to 400-500 seconds from the beginning of the experiment. 


3.4 Comparison with measurements 


We now compare the observed Euglena cell velocity histograms with the predicted 
shape of the distribution. The histograms along with the theory curves are pre¬ 
sented in Fig. All the sets could be fitted equally well. We do not show all of 
them as they practically overlap. The correlation coefficient is in the range from 
0.944 to 0.955 in different sets. We see that the cell motion matches very well the 
predicted shape for systems with active noise of high amplitude as shown in Fig. 

It has a characteristic strong peak around zero velocity. The fit is good in the 
whole velocity range. Moreover, we observe that the distribution changes with time 
during the experiment. The peak close to zero velocity becomes even stronger with 
time while the wings of the distribution, corresponding to the normally moving 
cells, get somewhat weaker, which is consistent with the increase of the fluctuation 
amplitude Dv (see curves with Dv = 0.3 and 1.0 in Fig. right panel). 

We attempted to interpret the latter observation using the depot model with 
active noise. The parameters of the motion are as shown in Table As we see 
from the data, the ratio of energy influx to the fluctuation amplitude q/{2Dv) 
remains nearly constant during the experiment. The ratio 'y/{2Dv) is decreasing 
by about 50%. Although we cannot resolve Dv independently, we can assume that 
the viscous friction coefficient 7 stays constant as it reflects just the shape of the 
cells. We determined it in our previous work for Euglena cells 7 = 1.2 x 10“^ 
N/{m/s) [36]. Then, the decrease of the ratio 'y/{2Dv) indicates the increase of the 
fluctuation amplitude by about 50%. At the same time, as the ratio q/(2Dv) does 
not vary much (the variation is 5% - 7%), corresponding to values of g = 1.8x 10“^® 
W to 2.7 X 10“^® W, we can infer the increase of the incoming power by ca. 45%. 
The ratio c/d is increasing most significantly during the experiment, the total 
growth being more than 100%. The propulsion velocity that is dehned in the 
depot model [2511111 by Eq. ( |19[ ) stays fairly constant at first, vg ~ 26 /rm/s, but 
decreases to 18 fxm/s towards the end. In the last interval the difference q/j — cjd 
turns negative and the propulsion velocity vg as parameter looses its meaning. The 
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Time label 

q/{2Dv) 

10-97/(2T)„) 

10®c/d 

Vq, fim/s 

100 

11.2 

5.6 

1.3 

25.9 

200 

10.7 

5.2 

1.3 

26.6 

300 

10.7 

3.8 

2.1 

26.9 

400 

10.7 

4.0 

2.2 

22.0 

500 

10.5 

3.3 

2.9 

17.7 

600 

10.4 

3.8 

2.9 

N/A 


Table 1 Parameters of the depot model with active fluctuations from fit of Eq. (l |24[ |) to 
Euglena gracilis one-dimensional velocity histograms. There is no detectable persistent motion 
of the cells for the given damping 7 at time 500 to 600 seconds as q /7 — c/d < 0. 


second term c/d expresses the rate of internal energy dissipation to energy transfer. 
Hence, the observation is consistent with the increase of internal dissipation in the 
cells during the experiment. The overall picture can be interpreted as follows: 
The Euglena cells exposed to intense illumination show negative phototaxis and 
are getting more and more fatigued at long exposures: The amplitude of velocity 
fluctuations as well are internal dissipation are growing while the fraction of energy 
that is converted into persistent motion is decreasing. 


3.5 Active multiplicative noise 

Here, we also demonstrate the application of the active noise concept to a model 
with multiplicative noise, which has been reported for random migration of the 
social amoeba Dictyostelium discoideum [8]. The authors of this study measured 
the first two moments of the velocity increments per unit time as function of the 
velocity parallel and perpendicular to the direction of motion. From those data 
they approximated the propulsion function and the noise term of a corresponding 
Langevin equation. 

In detail, they proposed based on their measurements, a simple linear friction 
term combined with a multiplicative noise term, with linear dependence on speed. 
Thus, the motion of the cells in their model is entirely driven by fluctuations with¬ 
out a deterministic component in the driving force. Please note that the authors 
considered in their model only passive noise, independent on the heading of the 
cell. 

Here we go beyond the friction model suggested in [8] and calculate the speed 
distributions of self-propelled agents with active multiplicative noise and so-called 
Schienbein-Gruler friction function [251117]. The function includes a constant driv¬ 
ing force along the heading direction 

ap{t) ='rvoeh(t). (25) 

and a linear damping by a Stokes friction —"/v (c.f. Eq. [^. Thus for no = 0 the 
model reduces to the situation observed in Dictyostelium migration [8]. However, 
for vq > 0, it also accounts for the case of individual agents with a finite preferred 
speed. Fluctuations along the heading direction depend on the velocity, whereas 
angular fluctuations are assumed to have constant intensity, i.e. 


v • 


’n^{t) = V2Dyg{v)^y{t)ei, + ^j2D^^^{t)e, 


(26) 
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Fig. 7 Stationary speed distribution p{s) Eq. \S0\ for the model with multiplicative noise 
Eq. \27\ for different values of Dy (left), 71 (center) and vq (right). Default parameters if not 
varied: 7 = 1 . 0 , 71 = 1 . 0 , Dv = 1 . 0 , r?o = l-O- 


Here, ^i{t) with i,j = v,ip are again sources of independent Gaussian white noise 
defined in Eq.([^. In agreement with [8] we consider the case 

= (1 +7l|?;|). (27) 


and V is the velocity component along the heading direction. Thus, the multiplica¬ 
tive noise grows along the heading direction in both directions. We would like to 
underline that the additive part in the noise for ?; = 0 together with the Gaus¬ 
sian character of this noise can perform transitions from a motion with a positive 
velocity along the heading direction to a motion with negative one and vice versa. 

Using the Stratonovich interpretation |48| we can now calculate the stationary 
speed distribution. The distribution pa{v,ip,t\vo,(fio,to) of the projected velocities 
V along the heading direction and on its angle ip obeys the Fokker-Planck equation 

+ + (28) 

The asymptotic stationary distribution \imtg^-caPa{v,p,t\vo,po,to) = pa{v) be¬ 
comes independent on the initial conditions and the angle. With continuity of the 
velocity distribution function density pa{v) at the origin v = 0 we obtain 


(29) 


Due to the lack of symmetry of the propulsion vq > 0, the probability density 
function in Pa('f) is also asymmetric with respect to the origin d = 0. Since the 
noise depends on the absolute value |i!|, it is symmetric and grows in negative and 
as well positive direction. This dependence induces different slopes for negative 
and positive velocities around v = 0 and pa is not differentiable at the origin for 
Vo > 0 . 

Summation of probability at positive and negative velocities v values yields the 
speed distribution 


p(s) oc (1 
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velocity Vx velocity Vx 



Fig. 8 Stationary velocity distribution of t he Ujj-velocity component P{vx,0), Eq. j3l| for 
the model with multiplicative noise Eq. l |27| for different values of Dv (left), 71 (center) and 
vq (right). Default parameters if not varied: 7 = 1.0, 71 = 1.0, = 1.0, vq = 1.0. 


Interestingly, the asymptotic behavior for large speed valnes decays exponentially 
and not as a Ganssian. In Fig. we show examples of the speed distribution given 
in Eq. (301 for different parameter sets. 

Finally, we transform the above results to Cartesian coordinates and obtain 
the two-dimensional Cartesian velocity probability density, which diverges at the 
origin: 


P(Va:,Vy)= ^ p[s=Jvl+vl) (31) 

The corresponding graphs of this function for Vy = 0 are shown in Fig. for 
different values of Dv, 71 and vq (as in Fig. [^. We stress that as in the case 
of additive noise, we obtain a pronounced peak of the velocity distribution at 
vanishing velocities. 


4 Conclusions 

In this work, we have discussed the concept of self-propelled particles with active 
fluctuations in the context of experimentally measured velocity distributions of 
Euglena gracilis algae. In biological agents, the locomotion is controlled by multi¬ 
ple internal decision processes, which cannot be easily accessed. These complex, 
unresolved influences will strongly contribute to the apparent stochasticity of the 
motion, and in the hrst approximation they can be conveniently summarized as 
active fluctuations. The intensity of such active fluctuations and its dependence 
on various environmental parameters yield important information on the inter¬ 
nal states of the biological agents and determines ecologically relevant quantities 
as for example the spatial dispersal (diffusion) of individuals |49| . Here, our re¬ 
sults show that the depot model [H] of active Brownian motion together with 
active additive fluctuations perfectly describes experimentally measured velocity 
distributions. The fitting of theoretical distributions to experimental data reveals 
a systematic change of locomotion parameters of individual algae with the dura¬ 
tion of the experiment, which provides important insights into the adaptation of 
individual cells to a stressful environment. 

We should add that the Euglena single-cell algae might demonstrate a more 
complicated behavior as it was reported previously. It might be one of the simplest 
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systems showing fatigue in response to external stress. Although we successfully 
described the population statistics, on the single cell level the mechanism of the 
active fluctuations and of the observed “fatigue” remains unclear. This mechanism 
can be investigated further if one addresses complex and not fully understood 
impulse-response behaviors involving the photoactivated adenylyl cyclase [391150) 
and protein kinase m- 

The theory of active fluctuations predicts a divergence of two-dimensional 
Cartesian velocity probability densities irrespective on the propulsion details |29| . 
The corresponding diverging probability density manifests itself by a finite (local) 
maximum at the origin, as it is the case in our experimental data. Similar maxima 
in Cartesian distributions have been observed also in very different experimental 
systems such as migrating Dictyostelium discoideum amoeba [9] or self-propelled 
colloids [52]. Thus we expect that active fluctuations are ubiquitous in living and 
non-living active matter systems and therefore should be taken into account in 
interpretation of corresponding experimental results and their theoretical descrip¬ 
tion. 
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